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I. INTRODUCTION 

With the recent discovery of a Higgs-like particle with a mass of about 125 GeV by both the 
ATLAS [1] and CMS [2] collaborations, the focus has now turned to deciphering the properties 
of this particle and determining whether it is the Standard Model (SM) Higgs particle, or part 
of an extended Higgs sector. Analyses performed by ATLAS and CMS collaborations have 
shown that the couplings of the newly discovered particle are consistent with a SM-likc Higgs 
boson, within the accuracy of their measurements. In light of the present precision of the 
Higgs data, the LHC can claim to have discovered a SM-like Higgs boson. However, there 
is still plenty of room for deviations from SM behavior of 0(10%). A SM-like Higgs boson 
is easily achieved in an extended Higgs sector in the decoupling limit, where the lightest 
scalar is identified as the observed SM-like Higgs boson, and the heavier scalars are somewhat 
separated in mass (e.g. with a mass scale above 350 GeV [3, 4]). 

Before the mass of the Higgs boson was known, upper bounds on the SM Higgs mass were 
obtained by requiring that the running quartic coupling parameter avoid Landau poles (LPs), 
i.e., the coupling was required to remain finite up to a given energy scale A [5-7]. Lower bounds 
were obtained by requiring that the scalar potential remain stable during renormalization 
group (RG) evolution [8-13]. That is, the scalar potential is bounded from below at all scales 
between the electroweak scale and A. These bounds were contingent on the assumption that 
no new physics beyond the SM (BSM) enters between the electroweak scale and A. Turning 
around the argument, the existence of a LP or an instability of the scalar potential at some 
energy scale A suggests that new BSM physics must be present at or below A. 

After the discovery of the Higgs boson, previously obtained bounds were updated using two- 
loop renormalization group equations (RGEs) in Ref. [14] and three-loop RGEs by Ref. [15]. 
The most recent analysis of Ref. [16] has shown that the SM scalar potential becomes unstable 
at a value of A well below the Planck scale, if the Higgs boson mass is smaller than 129.6 ± 
1.5 GeV. 1 Taken at face value, these results would further imply that we live in a metastablc 
vacuum that will eventually (and catastrophically) decay via tunneling into the true vacuum. 
However, the lifetime of the metastable vacuum is many orders of magnitude larger than the 
age of the universe [16, 17]. On the other hand, if the electroweak vacuum is absolutely stable, 
1 The quoted uncertainty takes into account the parametric uncertainty of m t and a s and the effects of 
unknown higher order corrections [16]. 
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then the recent LHC discovery of a 125 GeV SM-like Higgs boson requires the existence of new 
BSM physics at an energy scale below a scale of A ~ 10 9,5 GeV, where there is an uncertainty 
of about 1 in the exponent due to parametric uncertainties of rn t , a s and the Higgs mass [16], 
in order to avoid the metastability of the SM vacuum. 

Although the prospect of existence of new BSM physics is exciting, there is no guarantee 
that the scale of the new physics is close to the clectroweak scale. Nevertheless, arguments 
motivated by naturalness of the electroweak symmetry breaking (EWSB) mechanism suggest 
that BSM physics should be present at or near the TeV scale (see e.g., Refs. [18, 19] for a 
review and a guide to the literature). Many models of new physics have been proposed to 
address the origin of ESWB, and many of these approaches possess extended Higgs sectors. 
However, in such models one must specify the BSM physics in order to study the behavior 
of running couplings between the electroweak scale and some very high energy scale A. At 
present, there is no direct experimental evidence that the origin of the EWSB scale is a 
consequence of naturalness. Adding additional Higgs multiplets at or near the TeV scale 
by themselves does not address the origin of EWSB. Indeed, one could argue that it makes 
matters worse by adding additional fine-tuning constraints. Nevertheless, in this paper we 
shall accept the fine-tunings required to sustain an extended Higgs sector near the TeV scale. 
After all, we know that multiple generations exists in the fermionic sector of the Standard 
Model. Thus, we should be prepared for the possibility that the scalar sector of the theory is 
also non-minimal. 

Here, we shall focus on the two-Higgs doublet model (2HDM), which was initially proposed 
by Lee in 1973 [20] (for a review, see e.g. Ref. [21]). It provides a richer Higgs particle 
spectrum, namely three neutral scalars and a charged pair. The 2HDM admits the possibility 
of CP-violation in the scalar potential, both explicit or spontaneous. In the limit of CP- 
conservation, two of the neutral scalars are CP-even, typically denoted by h and H, (where 
rrih. < mu) while the other neutral scalar is CP-odd, denoted by A. We shall consider a very 
general version of the 2HDM that is not inconsistent with present data. Such a model must 
possess a SM-like Higgs boson (within the accuracy of the present Higgs data). In addition, 
Higgs-mediated tree-level flavor changing neutral currents (FCNCs) must be either absent or 
highly suppressed. These conditions are achieved if the non-minimal Higgs states of the model 
have masses above about 350 GeV and if the Yukawa couplings are aligned in such a way that 
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the neutral Higgs couplings are diagonal in the mass-basis for the neutral Higgs bosons. The 
most general 2HDM parameter space allowed by the present data is somewhat larger than the 
one specified here. Nevertheless, the restricted parameter space outlined above is still quite 
general and incorporates the more constrained 2HDMs considered in the literature. 

The existence of additional scalar degrees of freedom in an extended Higgs sector provides 
an opportunity to cure the vacuum metastability problem of the SM Higgs boson. However, 
by demanding no Landau poles and requiring a stable scalar potential at all energy scales up 
to the Planck scale, one imposes strong constraints on the parameter space of the extended 
Higgs sector. Investigations of this type have been performed in extended Higgs sectors prior 
to the discovery of the Higgs boson in Refs. [22-28]. With the discovery and identification 
of a SM-like Higgs boson, the question of the validity of extended Higgs sectors up to the 
Planck scale has become more focused. A number of authors have considered the stability 
properties of extended Higgs sectors with additional singlet scalar fields [29-32] and 2HDMs 
with constrained scalar potentials [33-35]. 

In this paper, we examine the theoretical consistency of the most general 2HDM between 
the electroweak scale and the Planck scale, using the one-loop RGEs of the model to investigate 
the possible occurrence of Landau poles and instability of the scalar potential. We focus on 
the decoupling regime of the 2HDM where the 125 GeV Higgs boson is SM-like [36, 37], and 
assume Yukawa alignment in the flavor sector [38] to avoid Higgs-mediated tree-level FCNCs. 
Our aim is to exhibit the allowed regions of the 2HDM parameter space that are free from 
both Landau poles and vacuum instability below the Planck scale. In particular, a 2HDM 
that satisfies these constraints does not require further BSM physics to stabilize the theory. 

One of the distinguishing features of the most general 2HDM is the fact that the two 
scalar doublet, hypercharge-one fields are indistinguishable. One is always free to define 
new linear combinations of the scalar doublets that preserve the kinetic energy terms of 
the Lagrangian. A specific choice for the scalar fields is called a basis, and any physical 
prediction of the theory must be basis independent. In our analysis, we employ a basis- 
independent formalism introduced in Ref. [39]. We consider the most general 2HDM scalar 
potential (which is potentially CP-violating) and the most general Yukawa sector, which 
introduces three additional independent 3x3 matrix Yukawa couplings. Without additional 
assumptions, the latter yields Higgs-mediated tree-level FCNCs, in conflict with observed 
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data. In order to circumvent this, we impose a “flavor alignment ansatz”, introduced in Ref. 
[38], which postulates that the independent matrix Yukawa couplings are proportional to the 
corresponding quark and charged lepton mass matrices. In this case one finds that, in the 
mass basis for the quarks and leptons, the matrix Yukawa couplings are flavor diagonal, and 
the Higgs-mediated tree-level FCNCs are absent. One way to achieve alignment in the Yukawa 
sector is to introduce a set of discrete symmetries which constrain the Higgs scalar potential 
and Yukawa couplings. The so-called Type-I and II 2HDMs [42], and the related Type X and 
Type Y 2HDMs [40, 41] provide examples of this type. Indeed, Ref. [43] showed that the 
flavor alignment is preserved under RGE running if and only if such discrete symmetries are 
present. The flavor alignment ansatz is more general, but requires fine-tuning in the absence 
of an underlying symmetry. 

This paper is organized as follows: In section II, we review the basis-independent formalism 
as applied to the 2HDM. In section III, we describe the Yukawa sector and present the flavor 
alignment model used in this analysis. In section IV, we present our numerical analysis of the 
mass bounds governing the lightest scalar, which are derived by requiring the stability of the 
2HDM potential and the absence of Landau poles in the scalar quartic couplings below the 
Planck scale. Our analysis employs both the one-loop RG running of the quartic couplings, 
along with an estimate of the effects of the two-loop corrections. In section VI, we present 
our conclusions. The one-loop basis independent RGEs are presented in Appendix A, and 
stability conditions on the basis-independent 2HDM potential are derived in Appendix B. 


II. BASIS-INDEPENDENT TREATMENT OF THE 2HDM 

II. 1. The Higgs Basis 

In a generic basis, the most general renormalizable SU(3)cxSU(2)iXU(l)y gauge-invariant 
2HDM scalar potential is given by 

V = + m\ 2 ($ 2 ^ 2 ) - [ml 2 ${$2 + h.c.] 

+ |Ai($I$i) + |A2( < f ) 2 < f > 2) + (^2^2) + A4(*^>{^>2) ($2^1) 

+ {1^5 ( < ^ > l^ > 2 ) + + -V( < £ > 2 < £ > 2)] ( < f > { < L 2 ) + h.C.} , 


( 1 ) 
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where $ 1 , $2 are two hypercharge-one complex scalar doublets. The two doublets separately 
acquire vacuum expectation values (vevs) (<h°) = v\/y/2 and ($ 2 ) — ^ 2 /a /2 with the constraint 
v 2 = |ui | 2 + \v 2 1 2 — (246 GeV) 2 . The parameters Ai., 2 , 3,4 and are real whereas As^y 

and m\ 2 are potentially complex. The 2HDM is CP-conserving if there exists a basis in which 
all of the parameters and the vacuum expectation values are simultaneously real. 

We shall adopt a basis-independent formalism as developed in Ref. [39], which provides 
basis-independent 2HDM potential parameters that are invariant under a global U(2) trans¬ 
formation of the two scalar doublet fields, <h a —» U a i$b (a, b = 1, 2). 

It is convenient to define the so-called Higgs basis of scalar doublet fields, 


Hi = 


H t \ _ u*$i +v *$ 2 

H°i 


Ho = 



~V 2 ®1 + Vi$2 


( 2 ) 


so that (H (*) = v/V2 and (H 2 ) = 0. The Higgs basis is uniquely defined up to a rephasing of 
the H 2 held, H 2 —> e lx H 2 . In the Higgs basis, the scalar potential takes the familiar form , 2 

V = Yi(H\Hi) + Y 2 {h\H 2 ) + [Y 3 H\H 2 + h.c.] + \Zi(H\Hi) 2 + \Z 2 {h\H 2 ) 2 


+z 3 (h\h 1 ) {h\h 2 ) + z a {h\h 2 ) {H\Hi) + {\Z,{H\H 2 ) 2 + [Z 6 (H\Hi) 


+Z 7 (h' 2 H 2 )] {h\H 2 ) + h.c.} , 


( 3 ) 


where Y A ,Y 2 , and Z 1 , 2 , 3,4 are real parameters and uniquely defined, whereas Y 3 and Z 5j6 j 
transform under a rephasing of H 2 , viz., [Y 3 , Z 6 , Z 7 \ —> e~ tx [Y 3 , z 6 , Z 7 \ and Z 5 —> e~ 2lx Z 5 . 
Minimizing the scalar potential then yields 


Yi = -lZiv 2 , Y 3 = -\Z 6 v 2 . 


( 4 ) 


The scalar potential is CP-violating if no choice of x can be found in which all Higgs basis 
scalar potential parameters are simultaneously real. 

The tree-level mass eigenstates of the neutral scalars can be obtained by diagonalizing the 

2 As discussed in Appendix A, the squared-mass and coupling coefficients, Yi, I 2 , and Zi <2 , 3,4 can be expressed 
as U(2)-invariant combinations of the scalar potential coefficients and the vevs, whereas Y 3 and Z 5 e 7 are 
U(2)-pseudoinvariant combinations of the scalar potential coefficients and the vevs that are rephased under 
a U(2) transformation [39]. 
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TABLE I: q ki as a function of the neutral Higgs mixing angles in the Higgs basis. 


k 

Qkl 

Qk2 

1 

COS 012 COS 013 

— sin 0i2 — i cos 0i2 sin 9 \ 3 

2 

sin 0i2 cos 0i3 

cos 0i2 — i sin 0i2 sin 0 i 3 

3 

sin 0i3 

i COS 013 


neutral scalar squared-mass matrix in the Higgs basis [44, 45], 

( 


\ 


Z\ Re(^e) -Im(Z 6 ) 

M. —v 2 \ Re(Z 6 ) l[Z 3 + Z 4 + Re(Z 5 )} + Y 2 /v 2 -|lm (Z 5 ) • (5) 

-Im(Z 6 ) -|lm (Z 5 ) \[Z 3 + Z 4 - R e(Z 5 )] + Y 2 /v 2 / 

The diagonalizing matrix is a real orthogonal 3x3 matrix that is parameterized by three 
mixing angles # 12 , # 13 , and 623 (details can be found in Ref. [45]). In terms of U(2)-invariant 
combinations of the mixing angles and scalar potential parameters, the squared-masses of the 
three neutral Higgs bosons, denoted by hi, h 2 and h 3 respectively, are given by [45], 

m k — | Qk 2 1 2 Y 2 + d 2 i. q^Zi + \\q k 2 \ 2 [Z 3 + Z 4 — Re(Z 5 e 2,6,a3 )] 


+Re(q k 2 )Re(q k 2 Z 5 e 2l<>23 ) + 2q kl Re(q k 2 Z (i e t023 )|, for k = 1,2,3, ( 6 ) 

where the q kl are invariant combinations of the mixing angles shown in Table I. It is convenient 
to choose a convention where m\ < m 2 < m 3 (which can always be arranged by an appropriate 
choice of neutral Higgs mixing angles). The squared-mass of the charged scalars is given by 


m 2 H ± = Y 2 + \Z 3 v 2 . 


( 7 ) 


II.2. Decoupling Limit 


The decoupling limit corresponds to taking the squared-mass parameter of the Higgs basis 
held H 2 large while holding the Higgs quartic coupling parameters fixed. In the perturbative 
regime, we take \Z % \ <0(1) and Y 2 S> v 2 . In this case [37, 45], 

,2 - 


sin 9 


12 


v 

% 


r^j 


sind 13 ~ O 


( 8 ) 




In addition, the decoupling limit requires that 


which implies that 


Im(Z 5 e- 2 "») ~ O (J-) , 


Re(Z 5 e- 2 "“) = -|Z 5 


( 9 ) 


( 10 ) 


The overall sign in eq. (10) [which is not determined by eq. (9)] is fixed in the convention 

where m 2 < m 3 . Using the above results in eq. ( 6 ) yields 

,2 - 


2 v 2 
m 1 = Zj\v 


1+0, y 


mi = Y 2 + \v 


2 \ 1 


1 .,2 


mn — Y 2 + 5 v 


Z 3 + Z 4 - \Z 5 + 0 — 


Z 3 + Z 4 + \Z$ 


Y 2 

v 2 

°'y. 


( 11 ) 

( 12 ) 

(13) 


At energy scales below Y 2 , the effective low-energy theory corresponds to the Standard Model 
with one Higgs doublet. Consequently, in the decoupling limit the properties of h\ approach 
those of the SM Higgs boson. The non-minimal Higgs states are roughly degenerate in mass, 
m 2 ~ m\ ~ m H± ~ U 2 , with squared-mass splittings of 0(v 2 ), 


mj-ml ~ \Z 5 \v 2 , 
m 3 - m 2 H ± ~ l(Z 4 + \Z$\)v 2 . 


(14) 

(15) 


In the decoupling limit of a general 2HDM, the tree-level CP-violating and flavor-changing 
neutral Higgs couplings of the SM-like Higgs state h\ are suppressed by factors of Oiy 2 /Y. 2 ). 
The corresponding interactions of the heavy neutral Higgs bosons (h 2 and h 3 ) and the charged 
Higgs bosons (hf ± ) can exhibit both CP-violating and flavor non-diagonal couplings. If Y 2 is 
sufficiently large, then FCNCs mediated by the lightest neutral scalar can be small enough to 
be consistent with experimental data. However, for values of Y 2 of order 1 TeV and below, 
tree-level Higgs-mediated FCNCs are problematical in the case of a generic Yukawa sector. 


III. YUKAWA SECTOR 


The most general 2HDM Yukawa sector, describing Higgs-fermion interactions, includes six 
Yukawa matrices (as compared to three in the SM). In a generic basis, the Yukawa Lagrangian 
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for the Higgs-quark interactions is given by eq. (Al). Following the discussion of Appendix A, 
we can re-express the Yukawa Lagrangian in terms of the quark mass-eigenstate fields [46], 

-C Y = U L {rf{*? + $&*) - D L tf{rf{*i + r%$ 2 )U n 

+U L K(n?*t + r^ ( l>t)D R + + r+ h.c., (16) 

where are 3x3 Yukawa coupling matrices and K is the CKM matrix. 

Using eq. (2), one can rewrite eq. (16) in terms of the Higgs basis scalar doublet fields, 

-C Y = U l (k u H? + p u H?)U R - D l K\k ii H f + p u H^)U R 

+U l K(k^H+ + p D] H+)D R + D l (k^H° + p d 'H°)D r + h.c., (17) 


where 3 


Q_viv? + v Irg 

K* = - 


Q i Q 
p Q = Z 1 ’ 2 '^ + vvr} i 


( 18 ) 

v v 

Note that p Q —» e~ lx p Q with respect to the rephasing H 2 —» e zx H 2 . Since (H®) = v/y/2 and 
{H®) = 0, it follows that the k u,d are proportional to the diagonal quark mass matrices, Mu 
and Mu, whose matrix elements are real and non-negative, 


M v = 


VK 


U 


V2 


= diag (m u ,m c ,m t ), 


M d = 


VK 


D 


V2 


= diag (m d ,m s ,m b ). 


(19) 


The Yukawa couplings of the Higgs doublets to the leptons can be similarly treated by replacing 
U —>• N, D —>• E, Mu —> 0, Mu —> Me and K —> 1, where N = (p e , ;/ ;i , u T ), E = (e, p, r) and 
Me is the diagonal charged lepton mass matrix. 

Since the Yukawa matrices p UJ: - E are independent complex 3x3 matrices, it follows that the 
Yukawa Lagrangian exhibited in eq. (17) generically exhibits tree-level Higgs mediated FCNCs. 
The off-diagonal elements of the p u ' D matrices are highly constrained by experimental data 
to be very small. As first shown by Glashow, Weinberg and Paschos (GWP) [47, 48], it is 
possible to naturally eliminate tree-level Higgs mediated FCNCs if, for some choice of basis 
of the scalar fields, at most one Higgs multiplet is responsible for providing mass for quarks 
or leptons of a given electric charge. In the 2HDM, the GWP condition is usually imposed in 
four different ways by employing the appropriate Z 2 discrete symmetry [40-42, 49, 50]: 


3 As noted in eq. (A6), the p® are U(2)-pseudoinvariant combinations of the Yukawa coupling matrices and 
the vevs, whereas the kP are U(2)-invariants. 
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1. Type-I Yukawa couplings: r/f = r/P = r/[ J = 0, 

2. Type-II Yukawa couplings: r/f = r/P = ryp = 0. 

3. Type-X Yukawa couplings: r/f = r/P = r/p = 0, 

4. Type-Y Yukawa couplings: r)Y = = 0. 

For example, it follows from eq. (18) that in the Type-I 2HDM, 

n U,D,L _ V ' U,D,L 

p ~v; K ' 

and in the Type-11 2HDM, 


( 20 ) 


y = Hi K u 
vi 


»; 


( 21 ) 


In light of eq. (19), the p F (F = U, D, L) are, in these cases, diagonal matrices in which case 
the neutral Higgs-fermion Yukawa interactions are flavor-diagonal at tree-level. 

If only phenomenological considerations are invoked in choosing the Higgs-fermion Yukawa 
couplings, then it is possible to consider the more general case of the flavor-aligned 2HDM 
introduced in Ref. [38]. In this model applied to the Higgs basis, one imposes the following 
conditions 


p u = a U K U , 


p D = a u k u , and p L = a L n L , 


n ,,d 


l.,l 


( 22 ) 


which generalize the Type-I and II results exhibited in eqs. (20) and (21). In eq. (22), the 
alignment parameters, a u,D,L , are arbitrary complex constants. 4 The flavor alignment con¬ 
dition shown in eq. (22) is not imposed by any symmetry, and is strictly unnatural (i.e., it 
can be achieved only by a fine-tuning of the model parameters). Equivalently, as observed in 
Ref. [43], the flavor alignment is preserved under RGE running only in the case of Type I, 
II, X and Y Yukawa couplings. Nevertheless, one can imagine the possibility of new dynam¬ 
ics above the electroweak scale that could be responsible for an approximately flavor-aligned 
2HDM. Thus, in our analysis we shall employ the more general eq. (22), which is sufficient for 
satisfying the phenomenological FCNC constraints. 5 

4 In practice, if the magnitude of the alignment constants are too large, then some of the Higgs-fermion 
Yukawa couplings will develop Landau poles below the Planck scale. In our analysis, we will determine the 


allowed regions of the flavor-aligned 2HDM parameter space where such Landau poles are absent. 

5 By choosing to work in the decoupling limit where Y 2 > 500 GeV, we ensure that FCNCs generated by 

one-loop radiative effects are not too large to be in conflict with experimental data (see e.g. Ref. [51]). 
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IV. RG STABILITY AND PERTURBATIVITY OF THE 2HDM 

Let us assume that the observed SM-like Higgs boson (with rrih ~ 125 GeV) is part of 
a 2HDM in the decoupling limit with a flavor-aligned Yukawa sector, with no other new 
physics present beyond the 2HDM below the Planck scale . 6 We shall examine whether there 
are regions of the 2HDM parameter space that yield a consistent model under RG running 
from the electroweak to the Planck scale. In general, two potential problems can arise in 
the RG evolution. First, Landau poles could arise from the divergence of the 2HDM qnartic 
scalar couplings and/or Yukawa couplings. Second, the 2HDM scalar potential conld become 
unstable at a higher energy scale. The case of Landau poles is fairly straightforward, although 
the precise energy scale at which they arise cannot be strictly determined, since it lies outside 
the perturbative regime of the RGEs. In practice, we shall consider that a Landau pole occurs 
when the relevant coupling exceeds 100 for some energy scale A < Mpi. Indeed, once such 
a large coupling is reached, it will very quickly diverge at an energy scale very close to A. 
In our analysis, we employ the one-loop RGEs for the quartic scalar couplings of the 2HDM 
in the Higgs basis given in Appendix A. These equations are strongly coupled, and thus a 
divergence in one quartic scalar coupling will cause a divergence in the rest. The leading 
effects of two-loop running will be assessed at the end of this section. 

In the SM, the requirement that the scalar potential is stable at all energy scales below the 
scale A is easily implemented. It is sufficient to require that the SM quartic scalar coupling is 
positive, i.e. Asm (A) > 0 for A > v. Requiring that the 2HDM scalar potential is stable at all 
energy scales below the scale A leads to a more complicated set of conditions. In the 2HDM 
with an unbroken, or softly broken, Z 2 discrete symmetry that sets A@ = A 7 = 0 in eq. (1), 
the stability conditions were first obtained in Ref. [53], 

Ai>0, (23) 

A 2 > 0 , (24) 

A 3 > — V AiA 2 , (25) 

A 3 + A 4 — | A, 5 1 > — \J Ai A 2 . (26) 

6 Incorporating light neutrino masses via the seesaw mechanism [52] with the mass scale of the right-handed 
neutrino sector assumed to be of order a typical grand unified scale has a very minor impact on the consid¬ 
erations in this paper. 
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However, in the case of a completely general scalar potential, the corresponding stability 
conditions are far more complicated (with no simple analytic form). Ref. [54] provides an 
algorithm for deriving the stability conditions for a general 2HDM, with no symmetry or CP 
assumptions imposed on the 2HDM scalar potential. In terms of the Higgs basis parameters, 
this algorithm is summarized in Appendix B. Except for special cases for the quartic scalar 
couplings, the corresponding stability conditions must be determined numerically. 

We now describe in detail the procedure used in our analysis. We assume that we are in the 
decoupling regime of the 2HDM, where the mass scale of the heavy Higgs sector is of O(Ah)- 
In light of eqs. (7), (12) and (13), we henceforth set A 2 H = Y- 2 . 

1. Start with the SM Higgs potential dehned at the scale of the 125 GeV Higgs boson. 

2. Use SM RG evolution to run the Higgs-self coupling parameter A and the fermion mass 
matrices up to the scale A H .‘ 

3. Match the one-doublet Higgs potential with the 2HDM potential by taking Z\ = A (Ah) 
and k f = \/2Mf(Ah)/ v (for F = U,D). This establishes the low energy boundary 
conditions. The effects of the lepton masses are negligible and have been ignored. 

4. Scan over all other 2HDM quartic scalar coupling parameters Z, and Yukawa alignment 
parameters a F (F = U,D). The latter fix the values of the p F (A h)- 

5. Run the 2HDM RGEs for the Z tJ n F and p F up to higher energies A. Check for stability 
of the potential at the scale A using the procedure summarized in Appendix B. 

6. Stop the running if a Landau pole is encountered or if the stability conditions cannot 
be satisfied. 

For the scalar sector, we scanned over the parameter space using 100,000 points, with 
\Zi\ < 0( 1), for i = 2,..., 7, to enforce the decoupling limit. These points were also subject 
to the constraint that they obey the stability conditions presented in Appendix B. Note that 
when \Zi\ <C 1 for i — 2,..., 7, we recover the SM Higgs sector. The choice of A h is subject 
' Starting the RG evolution at mz , we use a five flavor scheme to run up to to* and a six flavor scheme 
above mt . Running quark mass masses at mz and m* are obtained from the RunDec Mathematica software 
package [55], based on quark masses provided in Ref. [56]. For simplicity, the effects of the lepton masses 
are ignored, as these contribute very little to the running of the Zj . 
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to the condition A 2 H 3> v 2 , so that we are safely in the decoupling regime. Moreover, in order 
for the 2HDM to be distinguishable from the SM Higgs sector, A h should not be significantly 
larger than 0(1 TeV). We considered two different values, A h = 500 GeV and 1 TeV, although 
the allowed parameter regime in which the 2HDM remains consistent up to the Planck scale 
is not especially sensitive to the precise value of Ah in the desired mass range. In the case of 
A h = 500 GeV, it is plausible that the heavy Higgs boson states could be detected in high 
luminosity LHC running. Indeed, as we shall demonstrate later in this section, differences in 
the squared-masses of the heavy Higgs states can provide an important consistency check of 
this framework. 

The Yukawa couplings play a fundamental role in this analysis. As discussed in Section III, 
we have employed the flavor aligned 2HDM to describe the Yukawa sector, with random 
complex alignment parameters whose moduli were varied by several orders of magnitude. The 
evolution of the Yukawa couplings in the flavor-aligned 2HDM was first performed in Ref. [57]. 
Notice that the running of the Yukawa couplings can also generate Landau poles. Due to the 
large size of the top quark mass, at least one of the Yukawa couplings will be of order one 
at the electroweak scale, so that a Landau pole in the top-quark Yukawa coupling below the 
Planck scale can be generated by the RG running. The alignment parameters, unique for both 
the up and down quark sectors, were log random generated in such a way as to prevent such 
Landau poles in the running of the Yukawa couplings up to Planck scale. In the RG running, 
the initial value of the top Yukawa coupling was taken to be yt(‘m t ) = 0.94, corresponding to 
an MS top quark mass of m t (m t ) = 163.71 ± 0.9 GeV [56]. The non-occurrence of Landau 
poles then leads to the constraints 8 

10^1 <0.95 and \a D \ < 81.5, (27) 

as seen in Fig. 1. These results are quite consistent with those obtained in Ref. [57]. 

The effect of the alignment parameters in the one-loop quartic scalar coupling RGEs is 
to bolster the negative Yukawa terms, thereby further driving the quartic scalar couplings 
to be negative during RGE evolution. The influence of the Yukawa couplings in the scalar 
couplings RG evolution is dominated by yf terms (where y t is the top quark Yukawa coupling) 

8 For A h — 1 TeV, we find \a u \ < 0.97 and \a D \ < 84. The figure corresponding to Fig. 1 looks nearly 
identical, so we do not display it here. 
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FIG. 1: Distribution of absolute values of the flavor alignment parameters, for regions of 2HDM 
parameter space which remain valid up to the Planck scale, assuming A h = 500 GeV. 


in the one-loop /3-functions, where they provide a negative contribution. In this manner, the 
large size of the top quark Yukawa coupling tends to drive Z\ negative at large energy scales, 
thus provoking an instability in the potential. This will occur unless the starting point value 
(at the electroweak scale) of Z\ is large enough. Since Zi is directly related to the lightest 
CP-even mass in the decoupling regime, requiring the stability of the scalar potential between 
the electroweak scale and the Planck one therefore yields a lower bound on m^. Similarly, if 
the initial value of Z\ at the electroweak scale is too large, then a Landau pole will appear in 
the running of Z\ below the Planck scale due to the fact that the leading Zi contributions to 
the /3-functions of the quartic scalar couplings are positive, thereby driving the quartic scalar 
couplings to larger values as the energy scale increases. Preventing the occurrence of Landau 
poles thus establishes an upper bound on Z 1; and thus on rrih- 

Within the SM, these demands can only be satisfied up to the Planck scale by a rather 
narrow window of Higgs boson masses, which excludes the observed value of 125 GeV. As we 
shall now see, the complexity of the 2HDM scalar potential “opens up” that narrow window 
to include the known value of the Higgs mass. 
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FIG. 2: RG running of 2HDM quartic scalar couplings, with A h = 1 TeV. Red points correspond 
to parameter choices for which an instability occurs in the scalar potential; blue points indicate the 
presence of a Landau pole. The upper solid black line indicates the occurrence of a Landau pole in 
the SM. The lower solid black line indicates the limit for which the SM potential becomes unstable. 

V. NUMERICAL ANALYSIS 

V.l. Results from one-loop RG running 

Let us now compare the effect of the one-loop running of the SM scalar coupling, with 
its effect on the 2HDM quartic scalar couplings. The results of our calculations are shown 
in Fig. 2, which we now analyze in detail. The full 2HDM running begins at A# = 1 TeV, 
where the Z % for i = 2,... ,7 are chosen. 9 The red points in Fig. 2 correspond to choices 
of parameters Z t for which an instability of the potential occurred for a given higher scale 
A > A h- The blue points correspond to parameter choices for which a Landau pole occurred 
during the RG running at some scale A > A h- These results are to be compared with the 
corresponding results of the SM Higgs sector also shown in Fig. 2: the upper solid line indicates 
the maximally allowed value of rri^ to avoid a Landau pole and the lower solid line indicates 
the minimal value of rrih needed to avoid a negative SM quartic scalar coupling, at all energy 



9 The corresponding plot for A h = 500 GeV looks nearly identical, so we do not exhibit it here. 
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scales below A. We recover the well-known one-loop SM result that 140 < rrih < 175 GeV 
in order to preserve vacuum stability and avoid Landau poles in the running of the quartic 
scalar coupling at all energy scales up to Mpl [5-12] . 

The distribution of red and blue points in Fig. 2 has some interesting features. First, there 
are no blue points above the SM-Landau pole line. In fact, although the 2HDM scalar potential 
has several scalar couplings, their contributions to the 2HDM /3-functions are mostly positive. 
As such, when one of these couplings starts to become very large in its RG evolution, the others 
will not be able to counteract that growth, and a Landau pole is reached. Consequently, the 
upper limit for the quartic scalar coupling Z\ that controls the value of rrih hardly differs 
from the corresponding SM result. Second, note the appearance of many blue points below 
the SM-instability line. These correspond to Landau poles that occur for relatively low values 
of rrih, which is equivalent to low values of Z\. However, even though the initial value of Z\ 
at Ah may be small, the values of other Z { can be large, and thus Landau poles in these 
couplings can be generated, yielding those blue points below the SM instability line. 

The most interesting aspect of our results concerns the distribution of the red points, which 
correspond to the violation of one or more of the 2HDM stability conditions at the energy 
scale A. We see a great “density” of points around the SM-instability line. These points may 
be interpreted as regions of 2HDM parameter space that constitute small deviations from SM 
behavior. But the remarkable difference with the SM result is the appearance of many points 
below and to the right of the SM-instability line. For these points, the instability of the scalar 
potential occurs at a larger value of A for a given value of rrih as compared to the SM. Indeed, 
the full impact of the 2HDM on the RG evolution may be best appreciated by examining the 
rightmost boundary of Fig. 2 corresponding to A = Mpp. On this boundary, we find both blue 
and red points, for a range of Higgs masses from about 118 GeV up to 175 GeV. Thus we see 
that a range of 2HDM parameters exists for which it is possible to have a SM-like Higgs boson 
with a mass of 125 GeV, without that mass value implying an instability of the potential (or 
a Landau pole) between the electroweak and Planck scales. 

Let us now analyze more closely the region of parameter space for which the 2HDM is 
consistent up to the Planck scale. According to Fig. 2, only a narrow range of rrih (which 
corresponds to a narrow interval of values of Z\) is consistent with a 2HDM with a stable 
vacuum and no Landau poles from the electroweak to the Planck scale. Since the 2HDM 
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FIG. 3: Histograms of squared-mass differences of the heavy scalar states for Ah = 500 GeV. The 
left panel shows values of squared-mass difference between the two heavier neutral states. The right 
panel shows the values of the squared-mass difference between the lighter of the two heavy neutral 
states and the charged Higgs boson. The histograms correspond to 2HDM parameters for which 
there are no Landau poles and vacuum stability is satisfied at all energies below the Planck scale. 

quartic couplings are all coupled together in their RG running, it follows that the allowed 
ranges for all Z, } not only Z l} will likewise be quite narrow. This has interesting implications 
on the scalar mass spectrum. In fact, in light of eq. 14, the squared-mass splitting of the 
two heavy neutral Higgs states depends primarily on \Z 5 \. Likewise, eq. 15 shows that 
the squared-mass splitting of the heavier neutral Higgs boson and the charged Higgs boson 
primarily depends on Z 4 and \Z 5 \. Since the possible values of Z 4 and \Z 5 \ are restricted to a 
narrow range of values, it follows that the squared-mass splittings of the heavy Higgs states 
should also be strongly constrained. 

For a 125 GeV SM-like Higgs boson, we have evaluated the squared-mass splittings of the 
heavier Higgs bosons for 2HDM parameters that are consistent with a stable scalar potential 
and an absence of Landau poles up to the Planck scale. The histograms shown in Fig. 3 
exhibit the distributions in arbitrary units of the squared-mass difference between the two 
heavy neutral states (which is positive by definition) and the difference between the lighter 
of the two heavy neutral states and the charged Higgs pair, for A h = 500 GeV. Given the 
formulae in section II.2, all the heavy scalars have masses of order A h in the decoupling 















































18 



min 

max 

mean 

std. dev. 

(m| — ml)/v 2 

0.01 

0.26 

0.09 

0.05 

(ml - m 2 H± )/v 2 

-0.20 

0.11 

0 

0.05 

(m| — m 2 H± )/v 2 

-0.07 

0.19 

0.09 

0.04 


TABLE II: Squared mass splittings of the heavier Higgs bosons of the 2HDM with Ah = 500 GeV, 
for 124 < m/, < 126 GeV, for points that survive up to the Planck scale, using one-loop calculations. 



min 

max 

mean 

std. dev. 

(m| — ml)/v 2 

0 

0.29 

0.09 

0.05 

(ml - m 2 H± )/v 2 

-0.23 

0.12 

0 

0.06 

(mj - m 2 H± )/v 2 

-0.08 

0.19 

0.09 

0.04 


TABLE III: Squared mass splittings of the heavier Higgs bosons of the 2HDM with A h = 1 TeV, 
for 124 < m/,, < 126 GeV, for points that survive up to the Planck scale, using one-loop calculations. 

limit. The statistics of these histograms are summarized in Table II. If the 2HDM is valid up 
to the Planck scale, then the mass differences among the heavy Higgs states must be quite 
small. This presents a challenge for heavy Higgs searches at future colliders. It may be that 
such a spectrum could only be reliably determined at a multi-TeV lepton collider. Indeed, 
if the heavy Higgs spectrum could be determined at some future collider, it would provide a 
nontrivial check of the present framework in which the 2HDM is valid up to the Planck scale. 

The results shown in Table II are not particularly sensitive to the value of A h- For example, 
if A h = 1 TeV, then the distribution of possible squared-mass differences yields the results 
shown in Table III. Of course, in this case the corresponding mass differences are even smaller, 
and the separate discovery of each of these new scalar states at a future collider is even more 
challenging. 


V.2. The effects of two-loops RG running 


In the SM, the inclusion of the two-loop terms in the RGEs shifts the scalar potential in¬ 
stability boundary to a higher energy scale, which lowers the minimum Higgs boson mass that 
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FIG. 4: Left panel: Scale shift for converting the one-loop scalar potential instability boundary to 
the two-loop scalar potential instability boundary for a SM-like Higgs boson. Right panel: Higgs 
boson mass bounds in the flavor-aligned 2HDM, incorporating the scale shift shown in the left 
panel, assuming that A u = 1 TeV. Red points indicate an instability in the running; blue points 

indicate the presence of a Landau pole. 

is consistent with a stable scalar potential all the way up to the Planck scale. In particular, 
the results in Ref. [14] yield a minimal value of m/,. ~ 129 GeV for vacuum stability. More¬ 
over, given the currently observed value of 125 GeV for the Higgs boson, the SM vacuum is 
metastable under the assumption of no new physics beyond the Standard Model below about 
10 10 GeV. This means that the effect of including two-loop effects in the RG running lowers by 
about 10 GeV the minimal value of the Higgs mass that is consistent with vacuum stability. 

We expect that employing the full two-loop RG analysis for the 2HDM would provide a 
similar downward shift in the lower bound of Higgs masses that survive up to the Planck scale, 
as well as increase the fraction of points that survive. In practice, implementing this full two- 
loop procedure is computationally expensive. Instead, we present a procedure to estimate the 
two-loop RG results. Note that the stability curve for the SM scalar potential at two-loops 
is both shifted to a higher energy scale, and is less steep as a function of the Higgs mass, 
relative to the one-loop SM scalar potential stability curve. In essence, going from one-loop 
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min 

max 

mean 

std. dev. 

(m| — ml)/v 2 

0 

0.31 

0.11 

0.05 

(ml - m 2 H± )/v 2 

-0.23 

0.12 

0 

0.05 

(m| — m 2 H ±)/v 2 

-0.09 

0.23 

0.11 

0.04 


TABLE IV: Squared-mass splittings of the heavy Higgs bosons of the 2HDM with A h = 1 TeV, for 
rrih ~ 125 GeV, for points that survive to the Planck scale, using the two-loop extended procedure. 


to two-loops shifts the stability curve energy scale to a higher scale for a particular Higgs 
mass. From our one-loop SM calculations and the two-loop SM calculations of Ref. [14], we 
determine the energy scale shift of the SM scalar potential stability curves due to the inclusion 
of two-loop RG running. Taking A# = 1 TeV, the resulting scale shift function is shown in the 
left panel of Fig. 4, which then yields our “two-loop” result shown in the right panel, which is 
obtained by applying the scale shift to our one-loop calculation. This shift is applied only to 
those points in which the scalar potential became unstable, not for points that hit a Landau 
pole before the Planck scale. The upper bound on the SM Higgs mass due to the absence of 
Landau poles does not exhibit a similar shift from one-loop to two-loop calculations. As in 
the case of Fig. 2, the case of A h = 500 GeV yields nearly identical results. 

In our one-loop calculations, only 707 (or 0.707%) of the 100,000 points analyzed survive 
to the Planck scale in the 123 GeV to 128 GeV region. With the conversion shift and a 
double check that they satisfy the stability requirement for 2HDM quartic scalar coupling 
parameters, 1,371 more points reach the Planck scale for a total of 2,078 (or 2.078%) at the 
Planck scale, an increase of 94% relative to the one-loop results. With an increase in the 
number of points, the “two-loop” squared mass splittings of the heavier Higgs bosons for 
points that survive up to the Planck scale are given in Table IV. Comparing Tables III and 
IV, we see that there exists only slight differences in the squared-mass splittings of the heavier 
Higgs bosons when the approximate two-loop effects are included. Nonetheless, the increase 
in the number of points for which the model remains consistent all the way up to the Planck 
scale is according to what one should expect, in light of the observation that the two-loop 
contributions increase the stability of the SM potential. Thus, this quick estimate suggests 
that the 2HDM parameter space corresponding to a stable scalar potential and no Landau 
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poles in the RG running to the Planck scale is somewhat larger than the parameter regime 
identified in the one-loop analysis. In particular, given the observed Higgs mass of the 125 
GeV, there exists a robust region of the parameter space for which the validity of the 2HDM 
and the stability of the Higgs vacuum is preserved up to the Planck scale. 


VI. CONCLUSIONS 

The discovery of a SM-like Higgs boson with a mass rrih = 125 GeV has focused attention on 
the validity of the Standard Model at higher energies. Putting aside the question of the origin 
of the electroweak symmetry breaking (e.g., accepting the fine-tuning of parameters inherent 
in fixing the electroweak scale), one can ask whether the Standard Model is consistent all the 
way up to the Planck scale. Refined calculations of the radiatively-corrected scalar potential 
suggest that the Standard Model vacuum is at best metastable (and long-lived), with a deeper 
vacuum located at field values near 10 10 GeV, well below the Planck scale. 

Adding new degrees of freedom has the potential of ameliorating the problem of an unstable 
vacuum. In this paper we considered the two Higgs doublet extension of the Standard Model 
(2HDM) and examined the range of parameters for which the 2HDM is stable and perturbative 
at all energy scales below the Planck scale. Our aim was to make the minimal number of 
assumptions regarding the structure of the 2HDM required by the experimental data. Since 
the observed Higgs boson is SM-like (within the accuracy of the limited Higgs data set), 
we considered the 2HDM with the most general scalar potential in the decoupling regime. 
The Yukawa sector was treated using the flavor alignment ansatz, in which the second set of 
Yukawa matrices is proportional to the SM-like set at the electroweak scale to protect against 
tree-level Higgs-mediated FCNCs. Although the flavor alignment condition is not protected 
by a low-energy symmetry (except in special cases, which lead to 2HDMs of Types I, II, X 
or Y), it provides a more general framework which at present is consistent with experimental 
data. 

We scanned over the scalar potential parameters and the flavor alignment parameters to fix 
the boundary conditions at the scale of the heavy Higgs states. We then employed one-loop 
RGEs to run the 2HDM parameters up to the Planck scale, and required that no Landau poles 
are encountered, without generating an instability in the scalar potential. In contrast to the 
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Standard Model, it is possible to have a SM-like Higgs boson with a mass of 125 GeV while 
maintaining the validity of the 2HDM up to the Planck scale. We also presented a scheme to 
estimate the effects of the RG-rnnning at two-loops, by applying a scale shift seen in going 
from the one-loop SM scalar potential stability curve to the two-loop SM scalar potential 
stability curve. Such effects increase the number of points in the 2HDM parameter scan that 
survive Landau pole and stability requirements up to the Planck scale. 

The larger range of allowed values of rrih in the 2HDM (as compared with the SM) is a 
direct consequence of the fact that the 2HDM scalar potential contains more quartic scalar 
couplings than the SM, which increases the stability of the potential at all scales between 
the electroweak and the Planck scale. In contrast, we observed that the theoretical upper 
bound on in the 2HDM based on the non-existence of Landau poles up to Planck scale 
hardly differs from the corresponding SM behavior. This can be understood as follows. In 
the SM, the negative top Yukawa contribution in the quartic scalar coupling /3-function drives 
that coupling to negative values during RG running, unless its starting point is sufficiently 
large. In the 2HDM, even if the initial values of some of the quartic scalar couplings are 
small, and even though the top quark contributions to the /3-functions are still negative, other 
couplings are allowed to have large values, which (in some cases) counterbalance any putative 
instabilities arising due to RG running. The 2HDM scalar potential is thus comparatively 
more stable than that of the SM. 

Finally, we have obtained bounds on the square-mass differences of the heavier Higgs bosons 
in the parameter regime where the 2HDM remains valid up to the Planck scale. If the 2HDM 
is realized in nature, this could provide an important check of the consistency of the model. 
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Appendix A: One-Loop Renormalization Group Equations 


The one-loop RGEs for the SM used in this analysis are provided by Ref. [58]. The 2HDM 
one-loop RGEs in various bases are given in Refs. [43, 54, 59-61]. The one-loop RGEs found 
in the literature typically assume a 2HDM scalar potential with a Z 2 symmetry, (fq —> 

<f> 2 —> —4*2, to avoid FCNCs and/or are explicitly CP-conserving. Here, we present one-loop 
RGEs for the full 2HDM, using a basis-independent approach and making no CP assumptions. 

In a general 2HDM, the Higgs fermion interactions are governed by the following interaction 
Lagrangian: 


--2V = Q° Sy U% + q\ + e\ t.hf 0 )* E° r + h.c., 


(Al) 


summed over a, a = 1,2, where $i j2 are the Higgs doublets, = icr 2 $i, Q° L and E® are the 
weak isospin quark and lepton doublets, and U R , D ° R , E R are weak isospin quark and lepton 
singlets. [The right and left-handed fermion fields are defined as usual: ^r^l = Pr,l where 
P RL = 1(1 ± 75 ).] Here, Q° L , E®, U R , D R , E R denote the interaction basis states, which 
are vectors m the quark and lepton flavor spaces, and % ,r] 2 , 77, , rj 2 ,rj 1 , p 2 are 3x3 

matrices in quark and lepton flavor spaces. 

The neutral Higgs states acquire vacuum expectation values, 


($°> = ^ 


(A2) 


where v a % = 1 and v = 246 GeV. It is also convenient to define 


j 


(A3) 


where tyi = —621 = 1 and en = e 2 2 = 0. 

It is convenient to define invariant and pseudo-invariant matrix Yukawa couplings [39, 45], 


^ — V a'la > 


p f ’° - <€° , 


(A4) 


where F = U, D or E. Inverting these equations yields 

Va’° = K F ’°Va + p F ’°W a . (A5) 

Note that under the U(2) transformation, <f> a —» U a i&b [cf. eq. (A35)], 

n F ’° is invariant and p hfi —> (det U)p F '° . (A6) 
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The Higgs fields in the Higgs basis are defined by [45] 


H! = m a , 


H 2 = W*a$ a , 


(A7) 


which can be inverted to yield <f> a = Hiv a + H 2 w a . One can rewrite eq. (Al) in terms of the 
Higgs basis fields, 

--2V = W L + H 2 p u ’°) U° R + Q° l (. H lK D ’ ot + H lP D ’ ot ) D° r 

+E° l {H iK e ^ + H lP E ^) E° r + h.c., (A8) 

The next step is to identify the quark and lepton mass-eigenstates. This is accomplished 
by replacing Hi —> (0, v/y/2) and performing unitary transformations of the left and right- 
handed up and down quark and lepton multiplets such that the resulting quark and charged 
lepton mass matrices are diagonal with non-negative entries. In more detail, we define: 

PlU = V r P l U° , P R U = V E P R U°, P l D = V e P l D\ P r D = V e P r D\ 

P l E = V e P l E° , P r E = V R PrE° , P l N = V e P l N° , (A9) 

and the Cabibbo-Kobayashi-Maskawa (CKM) matrix is defined as K = V R V E ^. Note that 
for the neutrino fields, we are free to choose V R = V r E since neutrinos are exactly massless 
in this analysis. (Here we are ignoring the right-handed neutrino sector, which gives mass to 
neutrinos via the seesaw mechanism). 

In particular, the unitary matrices V[ and V R (for F = U, D and E) are chosen such that 


M v = -^=V£ r n Ufi V R ] = diag (m u , m c , m t ), 
M d = ~^=V R k d ’°^V r t = diag(m d , m s , m b ), 


M e = -^=V e k e '°W r ] = diag (m e , , m T ). 


It is convenient to define 


K c ' = YW B l ' , , k d = Vg K D ’°vF ', = V?k e ' 0 V?' 


R 


L * 


' R 


L > 


n u _ yUJJflyU] D _ yD n D,0yD f E _ yD EflyE] 

P ~ V L P V R i P — V RP V L i P ~ V RP v L 
Eq. (A6) implies that under the U(2) transformation, 4> a —> U ab & b , 

k f is invariant and p F —* (detf/)p F , 


(A10) 

(All) 

(A12) 

(A13) 

(A14) 

(A15) 
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for F = U, D and E. Indeed, n F is invariant since eqs. (A10)-(A12) imply that 


M F = — / 

V2 


(A16) 


which is a physical observable. The matrices p L \ p D and p E are independent pseudoinvariant 
complex 3x3 matrices. The Higgs-fermion interactions given in eq. (A8) can be rewritten in 
terms of the quark and lepton mass eigenstates, 


--2V = U l (k u H°' + p U H^)U R - D l K\k u H r + p u H 2 )U R 

+U l K{k d 'H+ + p d 'H+)D r + D l {k d 'H* + p d 'H%)D r 
+N l (k e *H+ + p E] Ht)E R + E l (k e 'H° + p E ^H° 2 )E R + h.c. (A17) 


We now write down the renormalization group equations (RGEs) for the Yukawa matrices 
Va ,0 i Va'° and V E ’°- Defining V = 16n 2 p(d/dp), the RGEs are given by [43]: 


= 


2 + b 2 + ^ /2 K’° + \ 3Tr[»7^°(#°) t + V^ty] + H 


,FO(E,o,n I c/,0 


-2(^’°)W + , (A18) 

^f’° = -h + y + y 2 w° +(3 Tr^f’W* 0 + } v r 


-2 + r,r(r>nw° + saw + °)V °, ( m 9 ) 

^f’° = -(^ 2 + ¥^ 2 )^’° + -°) V’ 0 + «’°) t e°] + Tr[(^°)t^] 

^ \vt°)V° ■ (A20) 


The RGEs above are true for any basis choice. Thus, they must also be true in the Higgs 
basis in which v = (1,0) and w = (0,1). In this case, we can simply choose r/, F '° = k, f, ° 
and r) E, ° = P F ’° 1° obtain the RGEs for the k f, ° and p F, °. Alternatively, we can multiply 
eqs. (A18)-(A20) first by v* and then by w*. Expanding %, which appears on the right- 
hand sides of eqs. (A18)-(A20), in terms of n' and p 1 using eq. (A5), we again obtain the 
RGEs for the cc F, o and p F, °. Of course, both methods must yield the same results, since 
the diagonalization matrices employed in eqs. (A10)-(A12) are defined as those that bring the 
mass matrices to their diagonal form at the electroweak scale. No scale dependence is assumed 
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in the diagonalization matrices, and as such they are not affected by the operators T>. 

Vk u ’° = -(8g 2 + f g 2 + g 9 ,2 )k u ’° + {3Tr [k u ’°k u ’ 0 ' + k d '°k d ’° t] + Tr [k e V^] |^° 

K u ’°p u ‘°t + « D ’°p D,0t ] + Tr[/^V’ ot ] jp 77 ’ 0 - 2(^’°t^’ 0 ^ 0 + p D ' ot ^V°) 
T/^’ 0 ^ 0 ^ 0 + p^tp^O) + i( K A 0 t K A 0 + p Ao tp Ao )/c r/,o + I(^, 0^,01 + p tf,o p iW )/c t/,o s 



Vp u ’° = — (8g 2 + fp 2 + gp'V’ 0 + "■V 7 ’ 0 * + p s ’°^’°t] + Tr [p 



^°^°tl UC/,0 


|3Tr [p^p^t + p D V iJt ] + Tr [p B V’ 0t ] \p U '° ~ 2(« D ’ ot p 


D,o K c/,o + p D,o tp D,o p c/,o) 


+p t/ ’°(^°t«^ 0 + p^tp^O) + + p D,o tp D,o )pC /,o + 1(^,0^,Of + p U,O p UM )p U,0 * 

Vk d » = -(8 g 2 + fp 2 + y 2 )n Dfi + |3TV[/c D - 0t /c D >° + k u ^k u ’°] + Tr [k e ’°'k e ’°\ }k d ’° 

+ <|3Tr [p D,ojl k d, ° + p u ’ 0 ' K u ’°] + Tr [p E,0 ^K E, °] j>p D, ° - 2(k d ’°k u ’°k u ^ + p D ’°K U ’°p u ^) 

+ (k D ’°K D ’ °t + p D ’V 0t )« Z? ’° + ±K D ’ 0 (K U ’°K Um + p U ’°p Um ) + Ik D ’ 0 (k D ’°'k D ’° + p D ’ 0t p D ’°) , 
Vp D, ° = -(8 p 2 + fp 2 + f^ ,2 )p D ’° + (3Tr[ K D ’°tp B ’° + K^tp^o] + 'n[« B >°tp®.0]\„A0 


+ |3Tr [p D ’°tp D -° + p^otpt/, 0 ] + Tr [ p E ’°t p E ’°]j> P D, ° - 2(« D ’ 0 p C7 ’ 0 K r7 ’°t + p D ’°p u ’°p u ^) 

+ ( K D ’°K D ’Of + p A0 p D,0 t)p D,0 + l p A0 (/c U ) 0 /c tt I 0t + p U,0 p tt,0 t) + l p A0 (/c A0 t/c A0 + pJ D,0 tpJ D, 0) ; 


= _Q p 2 + ^2^,0 + I 3Tr [^ot^,o + + Tr ^ot^ I e,o 


+ |3Tr [p D ^n Dfi + p u ’°U u ’°] + Tr [p E ’°\ E ’°] jp B ’° 

+(k e ’°k e ^ + p E ’V’° V’° + Ik E ’°(k e ^k e ’° + p E ^p E ’°), 


Vp E ’° = -(f g 2 + f g ,2 )p E ’° + | 3Tr [k d ^ P d ’° + k u ^p u ’°] + Tr^V’ 0 ] U E ’° 


+ |3Tr [p D ^p D ’° + p u ’°'p u ’°] + Tr [p B ’ 0 V’°] }p S ’° 

+{k e k e ^ + p E ’°p E ^)p E ’° + Ip e ’ 0 (k e ^k e ’ 0 + p E ^p E ’°) 
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Using eqs. (A13) and (A14), we immediately obtain the RGEs for the n F and p , 


Vk u = -(8 g 2 + | g 2 + g g' 2 )n u + |3Tr [ K U K U * + k d k d + Tr [k e k e ^] |k u 

+ |3Tr [k u p u ^ + n D p DJf ] + Tr [/sr^p^] |p^ — 2K(k 0 ^kP n u + p D ^ k d iPp u ) 

+ k U^ k U + p U\ p U ) + + p D\ p D )K \ K U + I (/ , V t + pUpU^ K U ; (A21) 

Dp^ = _(8p2 + |p 2 + gp' 2 ^ + |3Tr[pVt + + Tr[pV f ] jr/ 7 

+ |3Tr [pyt + p D p Dt ] + Tr[p E p Et ] }/^ - 2 X(/^tp^t^ + pP^ p ^K^p u ) 

+ pV^ + P^VO + \K{k d 'k d + p D 'p D )tfp u + \{k u k u t + p7 f )/ , (A22) 


I?^ = -(8p 2 + fp 2 + Ap' 2 )« D + 




+ +Tr ^t^] U^-2(^xWt + p ^t^ p t/t )A: 


+ ( k d k d t + p D p Dt )K D + + p u p u ^)K + \k d (k d ^k d + p Dt p D ), (A23) 


Vp D = -{%g 2 s + \g 2 + y 2 )p D + 


' D 'fP + K u 'p u ] +Tr [^t p ®] 


+ 3Tr [p D ^p D + p u ^p u ] + Tr [p*V] K p D _ 2 (« D ^pVt + frtp u p u ')K 


{k d k d ^ + p D p D] )p D + + p u p U] )K + ip D ( K Dt K D + p Dt p D ), (A24) 


^ £ = -(!p 2 + W 2 )^ 


k d\ k d + k u\ k u 


] +Tv[k e1 k e ]^k e + 


d ^k d 


+ P U] k u ]+Ty[p e \ e ] \p E + (k e k E] + p e p E] )k e + \k e (k E] k e + p E] p E ), (A25) 


?V ? = -(!p 2 + W 2 K + 


K Dt p D + + Tr[K Et p B ] \k E + \ 3Tr [p D] p D 


+ p u 'p u ] + Tr [p Et p E ] \p E + {k e k e t + pVV + + p £t p E ) • (A26) 



The 2HDM scalar potential in a generic basis shown in eq. (1) can be written in a more 
compact form following the notation of Ref. [39], 

V = Y a - b (&M + t $d) . (A27) 

Hermiticity requires that Y ab = Y t *- and Z abcd = Z* m -. In addition, the form of the scalar 
potential given in eq. (A27) implies that Z abcd = Z cdab . The full one-loop /^-function for Z abcd 
is given by, 

'D^'abcd ^^abef^cdfe T 2,Z a J c ^Zf bed 22/ a j ed Z j5 c g T ‘2Z abe jZ c gf ( j T 2,Z aa j b Z cde J 

-(3 g r2 + 9 g 2 )Z abcd + 1(3 g 4 - 2 g' 2 g 2 + g’ 4 ) 5 ab 5 cd + (3 g r2 g 2 )S ad 6 cb - 4N c Tr [rj%^g%f] 

+4(Tr [#r 7 ?] Z ebcd + Tr [rj^rjf] Z aScd + Tr [#//?] Z abed + Tr [gfgf] Z abcB ) . (A28) 

The squared-mass and coupling coefficients of the 2HDM scalar potential in the Higgs basis 
[cf. eq. (3)] can be written in the form of invariants or pseudoinvariants with respect to the 
U(2) transformations, <f> a —y U ab Q b , as shown in Ref. [39]. The three squared-mass parameters 
are given by 


n = UE®sei, r 2 = UiSjffi», n = u <a29) 

and seven coupling parameters are given by 

Zi = Z abcd v* a v b v* 5 v d , Z 2 = Z abc( j w* a w b w* 5 w d , (A30 ) 

Z?j = z abcd vl v b w* s w d , Z 4 = Z abcd w* a v b v* 5 w d , ( A31) 

Z 5 = Z abcd v a w b v* 5 w d , (A32) 

Zq = Z abcd v* a v b v* 5 w d , (A33 ) 

^7 = Z abcd v* a w b w* 5 w d . (A34) 

Note that under a U(2) transformation, v a —» U ab v b , whereas 

w a —> (det . (A35) 


Consequently, Y 4 , Y 2 , -Z’1,2,3,4 are real U(2)-invariants, whereas Y 3 , Z 5> qj are potentially complex 
U(2)-pseudoinvariants, which are rephased under a U(2) transformation, 

[Y 3 , z 6 , Z 7 ] (det Ur 1 [Y :il Z 6 , Z 7 \ and Z 5 -> (det U)~ 2 Z 5 . (A36) 
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Using the above results in eq. (??), the one-loop RGEs for the quartic scalar couplings in 

the Higgs basis are: 

VZ X = 12Z 2 + 4Z 3 2 + 4Z 3 Z 4 + 2Z 2 + 2jZ 5 j 2 + 24|Z 6 | 2 - (3 g' 2 + 9g 2 )Z 1 + |(</ 4 + 2 g’ 2 g 2 + 3 g 4 ) 

—4At.Tr [k Q] k q k Q] k q ] + 16(2Tr [k q U q ] Z 1 + Tr [k Q] p Q ] Z 6 + Tr [p Q ^ Q ] Z 6 *) , (A37) 

VZ 2 = 12Z 2 + 4Z 3 2 + 4Z 3 Z 4 + 2Z 2 + 2|Z 5 | 2 + 24|Z 7 | 2 - (3r / 2 + 9<? 2 )Z 2 + f(r / 4 + 2 g’ 2 g 2 + 3/) 

-4A c Tr [p Q] p Q p Q] p Q ] + 8 (2Tr [p Q] p Q ] Z 2 + Tr [k,^ p Q ] Z 7 + Tr [p Qt k q ] Z * 7 ) , (A38) 

PZ 3 = 2(Z 4 + Z 2 ) (3Z 3 + Z 4 ) + 4Z 3 2 + 2Z 2 + 2|Z 5 | 2 + 4|Z 6 | 2 + 4|Z 7 | 2 + 8 Z 6 Z 7 * + 8 Z 6 *Z 7 

- (3 g’ 2 + 9 g 2 ) Z 3 + ^ (g’ A - 2g r2 g 2 + 3<? 4 ) - 47V c Tr [k^kVV] + 4(2Tr [k q U q ] Z 3 
+2Tr [p Q] p Q ] Z 3 + Tr p Q ] Z 6 + Tr [p Qt k q } Z 6 * + Tr p Q ] Z 7 + Tr [p Qt k q } Z* 7 ) , (A39) 

25Z 4 = 2(Z 4 + Z 2 )Z 4 + 8 Z 3 Z 4 + 4Z 2 + 8 1 ^ 5 1 2 + 10|Zg|^ + 10|Z 7 | 2 + 2 ZgZ 7 + 2ZgZ 7 

-(3 g 2 + 9 g 2 )Z A + ^g ,2 g 2 - 4N c Trp Q p Q ^ k q ] + 4(2Tr [k q ^k q ]Z 4 + 2Tv[p^p Q ]Z 4 
+Tv[k q 'p q ]Z 6 + Tr[p Q ^K Q ]z; + Tv[^p Q ]Z 7 + Tr [ P q ^k q ]Z* 7 ) , (A40) 

VZ 3 = 2 Z 5 (Z 4 + Z 2 + 4Z 3 + 6 Z 4 ) + 10Z 6 2 + 10Z 2 + 4Z 6 Z 7 - (3 g' 2 + 9 g 2 )Z 5 - 4A c Tr [^p Q K^p Q ] 
+8 (Tr [k Q] k q ] Z 5 + Tr [p Q] p Q ] Z 5 + Tr [p Qt k q ] Z 6 + Tr [p Qt k q ] Z 7 ) , (A41) 

25Zg = 12Z 4 Zg + 6 Z 3 (Zg + Z 7 ) + 4Z 4 (2Zg + Z 7 ) + 2Zg(5Zg + Z 7 ) — (3 g 2 + 9g ,2 )Zg 
-4A c Tr [^k q k q ^p q ] + 4(3Tr Z 6 + Tr [p Q] p Q ] Z 6 + Tr [p Q ^n Q ] Z 4 

+Tr [p Q ^ Q ] Z 3 + Tr Z 4 + Tr [k q ^ p Q ] Z 5 ) , (A42) 

25Z 7 = 12Z 2 Z 7 + 6 Z 3 (Zg + Z 7 ) + 4Z 4 (Zg + 2Z 7 ) + 2Zs(Zg + 5Z 7 ) — (3 g' 2 + 9 g 2 )Z 7 


—4N c Tr [K Q ^p Q p Q] p Q ] + 4(3Tr [p Q] p Q ] Z 7 + Tr [k q ^k q ] Z 7 + Tr [p Q ^ Q ] Z 2 + Tr [p Q ^n Q ] Z 3 
+Tr [p Q ^ Q ] Z 4 + Tr [k Q] p Q ] Z 5 ) . (A43) 
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Finally, we note that the anomalous dimensions, which contribute to the quartic scalar 
coupling /3-functions, are given by 


'lab 


1 

32W 


(3 g 12 + 9 g 2 )S ab + 


1 

47r 2 




(A44) 


Appendix B: Bounded from below conditions for a general 2HDM potential 


To ensure the existence of a stable vacuum, the 2HDM scalar potential must be bounded 
from below, i.e. it must assume positive values for any direction for which the fields are 
tending to infinity. This places some restrictions on the allowed values of the quartic scalar 
couplings. For the case of the scalar potential given in eq. (1) with A 6 = A 7 = 0, those 
necessary and sufficient conditions are given in eqs. (23)-(26). 

We now review the analogous conditions for the most general renormalizable 2HDM po¬ 
tential, found in Refs. [54, 62]. It is particularly convenient to introduce a new notation for 
the scalar potential, based on gauge invariant field bilinears. Indeed, in many 2HDM studies, 
such as the comparison of the value of the potential in different vacua, the classification of 
scalar symmetries and stability conditions, the bilinear formalism provides a significant sim¬ 
plification in the calculations. This formalism also reveals a hidden Minkowski structure in 
the potential, which was established in Refs. [54, 62], A similar Minkowskian notation has 
been employed in Refs. [63, 64], 

There are four independent gauge-invariant field bilinears, which are defined by 


Tq — $1$! + < f > 2 < f > 2, 

n = - (<h][<f >2 + $ 2 $i) = -2 Re (V[$ 2 ) , 
r 2 = i - $ 2 $i) = — 2 Im , 

r 3 = - (V^i - $$$ 2 ) • 


(Bl) 


These four quantities form the components of a covariant four-vector, ry = (r 0 , r) with 
respect to SO(3,l) transformations. We also define A 1 = = (r 0 , —r) where g ,iu is the 

usual Minkowski metric. It is straightforward to verify that ry > 0 and Wry > 0, the latter 
being a consequence of the Schwarz inequality. That is, the four-vector r ;/ lives on or inside the 
forward lightcone LC + . The vacuum that preserves SU(2)xU(l) clectroweak symmetry [i.e., 
($ 1 ) = ($ 2 ) = 0] corresponds to the apex of LC + ] all neutral vacua correspond to the surface 
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of LC + , and any charge breaking vacua would lie on the interior of LC + . Transformations 
of the scalar fields that preserve the scalar held kinetic energy terms leave ro invariant and 
correspond to SO(3) rotations of the three-vectors, r. 

In terms of the bilinears defined in eq. (Bl), the scalar potential of eq. (3) can be written 
as 

V = -MX + XA/ry, (B2) 

with the 4-vector M /x and the mixed tensor A ^ given by 

M m = (-ffi + Y 2 ), Re y 3 , -Im Y :il ^-(Y, - y 2 )^ (B3) 

and 

\{%i + Z 2 ) + Z 3 — Re (Zq + Zj) Im (Ze + Zj) — |(Zi — Z 2 ) 

Re (Z 6 + Z 7 ) —Z± — Re Z 5 Im Z 5 —Re (Z 6 — Z 7 ) 

—Im (Z 6 + Z 7 ) Im Z 5 —+ Re Im (Z 6 — Z 7 ) 

UZ 1 -Z 2 ) -R e(Z 6 -Zr) Im (Z e - Z T ) -i(Z, + Z 2 ) + Z 3 

(B4) 

To ensure that the scalar potential is bounded from below one needs to evaluate the eigen¬ 
values and eigenvectors of the matrix AThen one can determine conditions on those 
eigenvalues and eigenvectors such that r^A^rv > 0. The eigenvalues A a (a = 0,1,2,3) of the 
matrix A ^ will be determined by the usual characteristic equation, 

det(A/-A a X) = 0. (B5) 

since g,^ = 5” is just the 4x4 identity matrix. The corresponding eigenvectors corresponding 
to eigenvalue A a will be denoted by V^ a \ For the most general 2HDM potential, the eigenvalues 
are the solutions of a quartic equation, which can in principle be determined analytically 
(although the corresponding expressions are not particularly transparent). However, it is 
straightforward to numerically evaluate the eigenvalues and corresponding eigenvectors. Note 
that, in general, some of the eigenvalues may be complex (since the real matrix A t / is not 
symmetric unless Z e = Z 7 = 0 and Z 7 = Z 2 ). 
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Having evaluated the eigenvalues and eigenvectors of A M ", we make use of Proposition 10 
of Ref. [54] to conclude that the 2HDM potential is bounded from below if and only if the 
following conditions are met: 

1. All the eigenvalues A a are real. 

2. Aq > 0. 

3. Ao > {Ai, A 2 , A 3 }. There may or may not be degeneracies among the three eigenvalues 
A* (* = 1,2,3). 

4. There exist four linearly independent eigenvectors corresponding to the four eigen¬ 
values A a , for a = 0,1, 2, 3. 

5. The eigenvector R(°) = (uooWioAboAbo); corresponding to the eigenvalue A 0 , is real and 
time-like. That is, it can be normalized so that 

| 4 0) | 2 = 4 ) - 4 - 4 - 4 ) = !• 


6 . The remaining three eigenvectors V W = (v 0 i, i’u, v 2 i, v 3i ) are real and space-like, i.e. 
normalized so that 


|yW | 2 = 


V 0 i - V U - v 2i - v 3i 


= - 1 . 


To illustrate this technique, we shall reproduce the bounded from below conditions for 
a potential with a Z 2 symmetry in the Higgs basis so that Zq = Z 7 = 0. Without loss of 
generality, we can choose Z 5 real by rephasing the Higgs basis field H 2 - The matrix A = A 
is then given by 


A — - 
2 


( \{Z 1 + z 2 ) + z 3 0 

0 — Z 2 — Z~, 

0 


\ 


(B 6 ) 


0 —\{Zi — z 2 

75 0 0 

0 — z± + z§ 0 

V K Z 1 - Z 2 ) 0 0 -\{Z l + Z 2 ) + Z 3 ) 

so that two of its eigenvalues can be immediately read off as Ai = — Z 4 — Z5 and A 2 = — Z4+Z5. 
The remaining two eigenvalues are 


A ± = Z 3 ± \JZ X Z 2 . 


(B7) 
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Since the eigenvalues must be real, if follows that 


ZiZ 2 > 0. 


(B8) 


A + is the largest eigenvalue and thus must corresponds to the time-like eigenvector. Hence, 
we identify A 0 = Z 3 + \J Z\Z 2 and A 3 = Z 3 — \JZ\Z 2 . Imposing the requirement that the 
scalar potential is bounded from below, it follows that the eigenvalues obtained above must 
all be real and obey the following inequalities: 


A 0 > 0 Z 3 > —\/Z\Z 2 

Ao > {Ax, A 2 , A 3 } Z 3 + Z 4 — |Z 5 | > —\JZ X Z 2 , 


(B9) 

(BIO) 


which are the Higgs basis equivalents of eqs. (25) and (26). The time-like eigenvector is 
pq°) _ ^ o, 0, y\ where the components x and y are related via the eigenvector equation by 


Z\ + Z 2 — \J Z\Z 2 

y = z, -z 2 

Since the time-like normalization condition implies that x 1 2 3 4 5 — y 2 = 1, we obtain 


(Bll) 


x 2 = 


(Z, - Z 2 ) 2 


4\/ Z\Z 2 [Z\ + Z 2 ) 

Thus we see that we must have Z\ + Z 2 > 0, which when combined with eq. (B8) yields 


(B12) 


Zi > 0 , z 2 > 0. 


Thus we recover the Higgs basis equivalents of eqs. (23) and (24). 
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